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ABSTRACT 

Pulsation frequencies of acoustic modes are calculated for realistic rotating stellar 
models using both a perturbative and a two-dimensional approach. A comparison be- 
t ween the two yields validity domains which are similar to those previously obtained 
in iReese. Lignieres &: Rieutordl (|2006l) for polytropic models. One can also construct 
validity domains based on polynomial fits to the frequencies from the two-dimensional 
approach, and these also turn out to be similar, thus further confirming the agree- 
ment between the perturbative and two- dimensional approach at low rotation rates. 
Furth ermore, as was previously shown in lEspinosa. Perez Hernandez fc Roca Cortes! 
(2004), adjacent frequencies in multiplets come close together, thus forming pairs. 
This phenomena, exclusive to two-dimensional calculations, is shown to be an unlikely 
explanation of the close frequency pairs observed in 5 Scuti stars. A systematic search 
for all close frequency pairs in the calculated spectrum is also carried out. The number 
of close frequency pairs is shown to agree with what is expected based on a Poisson 
distribution, but does not match the number or distribution of close pairs in stars like 
FG Vir. Furthermore, a lack of close frequency pairs appears at low rotation rates, 
where frequency multiplets do not overlap. Delta Scuti stars currently reported as 
having close frequency pairs do not fall in this interval. 
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1 INTRODUCTION 

The new space missions, CoRoT and Kepler, are substan- 
tially increasing the accuracy with which stellar pulsations 
are observed. A look at CoRoT's and Kepler's first aster- 
oseismic results shows the progress that has been made in 
the photometric detection of solar-like oscillations, as well 
as the substantial increase in the number of de tected modes 
in earlier type stars such as S Scuti stars |Michel et al.l 
|200S| , iKieldsen et~afl l201Ch . The high duty cycle and the 
long observational runs (up to 5 months for CoRoT and 
several years for Kepler) are key factors in this dramatic 
improvement. As a result, there is a drive on the theoretical 
side to achieve accurate computations of pulsation frequen- 
cies in stellar models. Although a fair amount of success 
has been achieved in slowly rotating stars, much work is 
still needed before being able to accurately model pulsation 
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modes in rapidly rotating stars. This is because at rapid 
rotation rates, the centrifugal force causes stellar deforma- 
tion thereby transforming what used to be a one-dimensional 
eigenvalue problem into one that is two-dimensional. This is 
all the more problematic as many stars rotate rapidly. 

In order to address the effects of rotation, two dif- 
ferent methods have been developed. The first is the 
perturbative method which starts from a non-rotating 
star and then includes various corrections which model 
the effects of rotation. These corrections are gener- 
ally expressed as a power series in terms of the rota- 
tion rate which is ass umed to b e a small parameter. 
Previo us wo rks include Saiol ( 198 J^) . iGoueh fc Thompson! 



l|l990l ) and iDziembowski fc Goodel i 19921) f o r the sec- 
ond order, and ISoufi. Goupil fc Dziembowskil (Il998l ) and 



iKarami et all (120051 ) for the third order. The advan- 
tage of this method is that it is computationally less 
expensive and mode classification is more straightfor- 
ward, thereby allowing an efficient study of an entire 
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spectrum of pulsation modes for a given model. How- 
ever this method is only valid for slow rotation rates. 
The other approach consists in solving directly the two- 
dimensional eigenvalu e system. Previ o us works which fo- 
cus on p -modes include Clementl dl98l[). Yoshida fc Eriguchl 
[|200lL lEspinosa. Perez Hernandez fc Roca Cortes! l2004). 
iLignieres. Rieutord fc Reesel |20od ') and Lovekin fe Deupred 
(2008). Although computationally more expensive, this 
method can handle pulsation modes at any rotation rate. 
A natural question is at what rotation rates does it becomes 
necessary to use a two-dimensional approach rather than a 
perturbative approach. 

A number of previous works have focused on answer- 
ing the above question in the cas e of acoustic modes. In 
iReese. Lignieres fc Rieutord] (|2006t ), validity domains were 
established for third order perturbative methods. The un- 
derlying stellar models were polytropic and the pertur- 
bative calculations were, in fact, pol ynomial fits to the 
non-p ert urbativjj] frequencies. Later on. lLovekin fc Deupred 
(2008) did comparisons using more realistic models. Once 
more the perturbative calculations were polynomial fits to 
the non-perturbative frequencies. Given the larger error bars 
on the calculations, perturbative calculations were consid- 
ered to be valid up to higher rotation rates. Other works 
include lOuazzanT et all (120091 1 where non-perturbative fre- 
quencies were compared with frequencies resulting, this 
time, from a perturbative analysis rather than a polyno- 
mial fit. Furthermore, the effects of avoided crossings were 
included in the perturbative calculations and shown to im- 
prove the_a^ree^r^nt_with two-dimensional calculations. Fi- 
nally, [Sudrez^et^l] £2010) investigated the effects of pertur- 
bative and non-perturbative frequencies of polytropic mod- 
els on echelle diagrams and various frequency separations. 
They showed that including avoided crossings in perturba- 
tive analysis becomes necessary for correctly predicting the 
behaviour of the small frequency separation, even at rela- 
tively low rotation rates. Nonetheless, a systematic study of 
the validity of perturbative methods for realistic models is 
still lacking. Given the accuracy of the new pulsation data 
coming from the C0R0T and Kepler missions, such a study 
is needed in order to correctly interpret this data. There- 
fore, in the present work, we will compare two-dimensional 
calculations with frequencies resulting from a second order 
perturbativ e analysis an d estab lish validity domains similar 
to those in IReese et al.l (|2006t ). These calculations will be 
carried out for acoustic modes, using realistic stellar models 
rather than polytropic ones. Furthermore, we will discuss 
some of the qualitative differences between the two types of 
calculations. 

The next section deals with the stellar models. This is 
then followed by a description of the perturbative and two- 
dimensional pulsation calculations. Section [4] compares the 
results from the two approaches and section [5] deals with 
close frequency pairs, both within multiplets, and for an 
entire spectrum. Finally, a discussion concludes the paper. 



1 The term "non-perturbative frequencies" in this context means 
frequencies based on a two-dimensional calculation, rather than 
frequencies in which the effects of rotation have been removed. 



2 MODELS 

Stellar models were c reated using the Aarhus stella r evo- 
lution code, ASTEC (IChristensen-Dalsgaardl l2008b[) . The 
OPAL 1995 opacity tables l|lglesias fc Rogers! I1996T ) wore 
used in co njunction with the Kurucz low temperature ad- 
justments (|Kuruczlll99ll ). and the Eggleton, Faulkner and 
Flann ery equation of state l|Eggleton. Faulkner fc Flannervl 
Il97&t ). Evolution begins at zero-age main sequence (ZAMS) 
and continues to a user-specified age. The initial mass, M, 
angular velocity f2 and hydrogen and heavy element abun- 
dances X and Z are also set by the user. The models cal- 
culated here are IMq and 2Mq ZAMS models with so- 
lar elemental abundances and rotation rates ranging from 
Q/D.K = - 0.25, where Q K = y/GM/R 3 and R is the stel- 
lar radius. A selection of the models calculated is given in 
Tabled] 

The effects of rotation were then included to second or- 
der in rotational velocity, Q,. In particular, the stellar struc- 
ture is no longer spherically symmetric but oblate due to the 
centrifugal force. We therefore perform a transformation of 
the coordinate system (r, 9, (j>) —¥ (x, 6, <f>), from a system of 
shells of constant radius, r, to shells of constant pressure, 
chosen such that x — R maps the surface of the star which 
can now be treated as a surface in hydrostatic equilibrium. 
The new radial coordinate x takes the form 



x = [l + hn(r)P 2 (cos8)]r, 



(1) 



where hn is the transformation coefficient and P2(cos#) = 
3cos 2 8-1 the I = 2 Legendre polynomial. This can be in- 
verted to give, to order fi 2 



r = [l-hn(x)P 2 (cose)]x. 



(2) 



We define a function u — hnx such that u satisfies the fol- 
lowing second order differential equation: 



JfTu = &(x), 
where 

X = 
<S{x) = 



d 2 . f8nx 2 p 2\ d i 



Ax 2 



+ 



Mq 



x I Ax x 2 ' 



4ttx po 



(;m^ \ fr2 - dl (xj ' - ! + 



6.r 



GMq 



(3) 

(4) 
(5) 



1 d / 2 d 
~G~MqoA~x \ X Ax- {xf62) 

and f r o = — |n 2 a;, f r 2 = |^ 2 a;, fe2 = \Q 2 x and go = 
m/M is the fractional mass interior to a radius x. Boundary 
conditions for this system are 



u regular as x — > 
U ' + -x=-Gii4x^ + ^ 



at x = R. 



The surface boundary condition is found by matching the 
gravitational potential <f>o onto the vacuum potential for x > 
7?. 

This change in coordinates leads to a change in all 
variables with a radial dependence such that e.g. po(x) = 
po(r) + pn(r)P2(cos9), where the subscript "0" represents 
the equilibrium quantity, with similar expressions for po(x) 
and cq (a;). 
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Table 1. Model parameters. 
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3 PULSATION MODE CALCULATIONS 
3.1 Perturbative analysis 

In order to calculat e perturbative frequencies we follow 
closely the method of lGough fc Thompson! (|l990l ). The per- 
turbation equation, to order fi 2 , can be written: 

j£f i + P quj 2 (; = wJt,t + Jf£ (6) 

where, 

set = -v [(p - poc 2 ) v ■ i- £■ vpo] + po v(v ■ i) 

V(lnp )Vp , 

= -2ip w-VC, (7) 
= [-£• V(u- W) + («. V) 2 ^ . 

Here po, po and Co are the equilibrium pressure, density 
and sound speed respectively, uu is the angular pulsation fre- 
quency and vo the velocity field resulting from rotation. £ is 
the Lagrangian displacement and is defined such that 

Wm = [£{r)Y e {e,(p),r){r) — , -— — (8) 



where £, r\ are the radial and horizontal amplitude functions 
respectively and Y/™ is the spherical harmonic for a mode of 
degree £, radial order n and azimuthal order m. We choose 
to normalise Y ; m such that 



|Y £ m (0,0)| 2 sin0d0d(£= 1, 



(9) 



and we normalise f(r) such that (,(R) = 1. 

The structural changes of the equilibrium model re- 
sulting from rotation cause perturbations to the eigenfunc- 
tion £ and to the function Jzf£, denoted respectively by 
£i and Ji?n£- For further details we refer the reader to 
iGough fc Thompson! (|l99Ct ). These perturbation quantities 
are substituted into equation (|6} and we linearise the re- 
sulting expression to second order in rotational velocity S7. 
This leads to the following expression for the oscillation fre- 
quency, valid to second order in f2: 

u n i m = u>o + i*i + (2/cjo) _1 {^o ■ (^6 — -^h — Pf2^o)£o) 
-cj 2 21 (2cj )~ 1 - wm/ _1 (/5oCo • £i) 



+ (2/)- 1 (£ -^bli) 

+wm(2Jwo)" 1 <fS -^ofo), (10) 

where (...) = j 2 ^ J fl ...x 2 sin 9dxd6d<f> is, to first order, 
the volume integral over the interior of the star. 6Jq is the 
unperturbed oscillation frequency, 

rR 

1= / p x 2 [C 2 + 1(1 + l)v 2 ] dx, (11) 

and u>o.i contains the terms to first order in f2 and can be 
expressed as: 



um = — flp x 2 V ) 2 + (£(£+ 1) - 2) V 2 ] dx. (12) 
1 Jo 

The ADIPLS adiaba tic pulsation package 
(jChristensen-Dalsgaardl l2008al ) was then modified to 
include second order rotational effects, and subsequently 
used to calculate perturbative frequencies. 



3.2 Two-dimensional calculations 

In order to carry out the two-dimension al pulsation c alcula - 
tions, we used the same approach as in iReese et al l l|2009h . 
The pulsation equations are expressed in terms of a new co- 
ordinate system which follows the shape of the star, then 
projected onto the spherical harmonic basis, and finally 
solved using the p ulsation code TOP (Two-dimensional Os- 
cillation Program. IReese et al.ll2009T ). 

The new coordinate system can be represented by 
(£, 8, <f>) where £ is the radial coordinate, 8 the colatitude, 
and (j> the longitude. The relationship between £ and r, the 
distance from the origin, is given by the following two equa- 
tions: 

a5 



r(C,0) = (l-e)C + 



5C - 3C J 



(R* s (e)-l + e) : 



(13) 



r(C,0) = 2e + (l-e)C 



+ (2C 3 - 9< 2 + 12< - 4) (Rt(8) - 1 - e) , (14) 



R 



R 



P ° 1 



where e = 1 - ^jf, R*(8) = ^ - P 2 (cost 
is the polar radius at the last grid point and R* the radial 
coordinate at the last grid point. The first equation applies 
to the first domain, i.e. for £ £ [0, 1], which corresponds to 
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the star. It must be noted that contours of constant x and 
£ values do not in general coincide, except at the last grid 
point where x = R* and £ = 1, because only the former 
corresponds to isobars. As a result, the model must be in- 
terpolated onto this new grid. The second equation applies 
to the second domain, i.e. for C, £ [1,2], which lies outside 
the star and in which only the perturbation to the gravita- 
tional potential is used. This coordinate system ensures that 
the standard regularity conditions can be used in the centre 
since it behaves like the spherical coordinate system in the 
centre, and it allows the use of simple boundary conditions 
both on the stellar surface for the Lagrangian displacement 
and pressure fluctuations, and on the outer boundary of the 
second domain for the perturbation to the gravitational po- 
tential. 

The pulsation equations are expressed in terms of the 
Lagrangian displacement, £: 

= p+V-(p f), (15) 

= [—uj + mQ] p £ — 2i [— lj + mQ] pofi x £ 

-Vp + pgeB-poVV, (16) 

= p + fV Po -c 2 (p + fVp o y (17) 

= A*-Ap, (18) 

where p, p and ^ are the density, pressure and gravitational 
potential fluctuations. This is th e same set o f equa tions as 
what is given in the appendix of I Reese et all l|2009l ) except 
that there is no gradient of the rotation rate, since the ro- 
tation profile is uniform, and each occurrence of uj has been 
replaced by — u, so as to ensure that prograde modes corre- 
spond to positive values of m, the azimuthal order. Further- 
more, in order to increase the accuracy of the calculations, 
the effective gravity has been calculated from the gradient of 
the total potential rather than from the gradient of the pres- 
sure divided by the density. For explicit expressions in terms 
of the spheroidal coordinate system, we refer the reader to 
iReese et all (|2009l ). 

Calculatin g pulsation modes in models based on the 
ASTEC code (|Christensen-Dalsgaardl l2008bh is similar to 
calculating modes in m odels based on the Self -Consistent 
Field (SCF) method |MacGregor et alj 120071 ) . although 
there are some noteworthy differences. ASTEC models in- 
clude atmospheres which means that the associated radial 
grids become very dense near the surface. As a result, care 
must be taken when expressing the radial differential op- 
erator in algebraic form. For insta nce, using the four th or- 
der finite difference implemented in Re ese et al. I ll2009h leads 
to numerically unstable behaviour. This has been replaced 
with a more stable version of fourth order finite differences, 
in which the differential equations are solved on carefully 
chosen intermediate grid points. 



4 VALIDITY DOMAINS FOR PERTURB ATIVE 
CALCULATIONS 

Before establishing validity domains for second order pertur- 
bative methods, it is important to compare the two meth- 
ods for non-rotating (or nearly non-rotating) stars. Differ- 
ences arise from the different formulations of the pulsation 




500 1000 1500 2000 

Frequency, / (in /iHz) 



Figure 1. A comparison of m = modes for a 1 Mq model 
rotating at Q = 10 — 5 S2k- The radial orders are n = 1, 10 and the 
harmonic degrees i = 0, 3. These differences are larger than 0.08 
/iHz, the error bar on a CoRoT long run, and must therefore be 
taken into account when constructing validity domains. 



equations and should be correctly characterised before pro- 
ceeding to establish validity domains. Figure [TJ shows the 
differences between m = modes calculated with the two 
methods. Although this error is small, it is larger than the 
error bar from a CoRoT long run and should therefore be 
taken into account when constructing validity domains. 

Figure [2] shows a comparison between the two methods 
for three different error bars. The first two correspond to 
a long and short CoRoT run and the third corresponds to 
2.3 days of observation. The radial orders of the modes are 
n — 1 — 10, the harmonic degrees 1 = — 3 and the azimuthal 
orders m — — I to I. Frequency multiplets calculated using 
perturbative theory were shifted using the differences which 
are plotted in Fig. [TJ This leads to somewhat larger and 
more realistic validity domains, especially for the smallest 
error bar. Based on these validity domains, non-perturbative 
effects will start to play an important role beyond « cq = 
30 km.s -1 for a long run, and v cq = 50 km.s^ 1 for a short 
run, in a 1 Mq star. However, given that solar-like pulsators 
tend to oscillate with radial orders between 15 and 25, these 
limits are expected to be lower, as based on the trends which 
can be seen in Fig. [5] 

If the perturbative frequencies are replaced with poly- 
nomial fits to the non-perturbative frequencies, similar va- 
lidity domains are obtained, as can be seen in Fig. [3] This 
comparison shows that there is a good agreement between 
the polynomial coefficients and the coefficients deduced from 
perturbative theory, thus providing a further check that both 
the perturbative and two-dimensional approach agree at low 
rotation rates. 

It is then interesting to compare these validity domains 
with those obtained from polytropic models. Figure [4] shows 
such a comparison. It includes the validity domain for N — 3 
and N = 1.5 polytropic models and compares them with 
that of the 1 Mq ZAMS model using the 0.6 pHz error 
bar. The polytropic models are scaled so as to have the 
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Figure 2. Validity domains for second order perturbative fre- 
quencies of 1 Mq ZAMS models for three different error bars, the 
first two correspondi ng to CoRoT error bars. These figures are 
analogous to Fig. 4 of|R eese et al 1 l|2006l l. which were established 
for an TV = 3 polytropic model. 



Figure 3. Same as Fig.[2]except that the perturbative frequencies 
have been replaced with second order polynomial fits to the non- 
perturbative frequencies. 

same mass and approximately the same radijf] as the 1 Mq 
ZAMS model. In order to superimpose the different valid- 



The polar and equatorial radii were made to satisfy "' p ° 1 ^"" Rcq = 
0.89 Rq, which is approximately true of the 1 Mq ZAMS model. 
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Figure 4. Comparison between between the validity domains 
for two polytropic models and the 1 Mq ZAMS model, for the 
0.6 /^Hz error bar. The threshold separating the regions where 
perturbative methods are valid and not valid is plotted so as to 
allow the superposition of the different domains. 
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ity domains, we have simply plotted the threshold between 
the region where the perturbative approach is valid, and the 
region where a two-dimensional calculation becomes neces- 
sary. As can be seen from the figure, the validity domains are 
quite similar. At high frequencies, they follow the same ten- 
dency as a function of frequency and have a similar spread 
around the mean value. At low frequencies, the threshold for 
the 1 Mq model seems to be, on average, lower than that of 
the polytropic models. This effect is, however, likely to be 
an artifact due to the cutoff at f2 = 0.25 Qk for the 1 Mq 
frequency calculations (see Fig. 



Figure 6. Upper panel: frequency differences between adjacent 
modes in the I = 1 — 3 frequency multiplets of 2 Mq models, 
as a function of the rotation rate. The frequency span of the 
spectrum is approximately 250 — 1060 /iHz. The grey band in the 
middle corresponds to differences which are less than 0.1 c/d. 
Lower panel: number of frequency differences below 0.1 c/d as a 
function of the rotation rate. This number is generally quite low 
except for slow rotation rates, where the average spacing between 
frequencies in a multiplet is comparable to 0.1 c/d. 



5 CLOSE FREQUENCY PAIRS 

In lEspinosa et all \2004 ), it was pointed out that adjacent 
frequencies within a multiplet tend to pair up and come 
close together, except for the most retrograde mode (i.e. 
I — —\m\ sectoral mode). It was then suggested that this 
p henomena could explain t he close fr equency pairs observed 
inlBreger fc Bis chof (200 ^) (see also [ Breger fc Pamvatnvkhl 
2006b and lBreger &: Pamvatnvk"hll2006a |). A si milar pairing 
up of modes also occurs for the frequencies in iReese et al.l 
(2006) as well as the non-perturbative frequencies presented 
here. In what follows, we will consider 2 Mq ZAMS models, 
as this is more representative of 8 Scuti stars. The second 
half of Table [T] gives the characteristics for a selection of 
these models. 

Figure [5] shows four sets of / = 2 multiplets calculated 
using both perturbative and two-dimensional calculations. 
As can be seen in the figures, only the two-dimensional cal- 
culations lead to this behaviour. The perturbative calcula- 
tions produce, instead, spacings which decrease uniformly 
when going from the most retrograde to the most prograde 
mode. However, the effects of avoided crossings have not 
been included in the perturbative calculations, so it remains 
to be seen whether including this effect can also produce a 
pairing up of adjacent frequencies. 



It is then interesting to investigate whether this ef- 
fect could explain the close frequency pairs observed in 
iBreger fc Pamvatnvkhl (|2006bh . In Fig. [6] we compare the 
frequency differences of mode pairs to a frequency separa- 
tion of 0.1 c/d (i.e. 1.1574 /iHz), a typical threshold for 
observed frequency pairs. As can be seen from the upper 
panel, these differences remain larger on average than 0.1 
c/d. The lower panel shows the number of sufficiently close 
frequencies as a function of the rotation rate. As can be 
seen, this number is not very large for any given rotation 
rate, except when the rotation rate is around or below 0.1 
c/d. However, 18 close frequency p airs were observed in FG 
Vir l|Breger fc Pamvatnvk h 2006a), and it's equatorial ve- 
locity is estimated to be 66 ± 16 k m.s , as based on the 
modelling of line profile variations l|Zima et al.| [2006). As a 
result, the pairing up of adjacent modes in frequency mul- 
tiplets do not seem to account for close frequency pairs in S 
Scuti stars. 

Of course, it is always possible to look at all of the fre- 
quency differ ences for a given spectrum. This app roach has 
been used bv lLenz. Pamvatnvkh fc Bregerl (|2008h to study 
3 S Scuti stars, including FG Vir. Here, we will push the 
analysis slightly further by looking at how the number of 
close frequency pairs depends on the rotation rate. Figure [7] 
shows the number of close frequency pairs as a function of 
the rotation rate, as well as where they occur. The dashed 
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Figure 5. Four £ = 2 multiplets calculated using perturbative (left) and two-dimensional (right) calculations in 2 Mq models. Frequencies 
only pair up in the two-dimensional calculations. In the perturbative case, the spacings decrease uniformly when going from the most 
retrograde to the most prograde mode. 



line in the lower panel gives the expected number of pairs, 



iV CX p., as based on a Poisson distribution: 



1) 



exp 



1)5/ 



A/ 



(19) 



where JV mo de is the number of modes in the spectrum, A/ 
the frequency span of the spectrum, and Sf the target fre- 
quency separation (i.e. 0.1 c/d). At a sufficient rotation 
rate, the average number of close frequency pairs matches 
A^cxp. indicating that the spectrum is behaving like a Poisson 
distrib utio n. As was pointed ou t bv | Lignieres fc Georgeotl 
l|200ctl and lLignieres fc Georgeotl I^OOgT the frequency spec- 
trum of rapidly rotating stars is subdivided into classes of 
regular modes, such as these low degree modes, and chaotic 
modes, the frequencies of which follow Poisson and Wigner 
distributions respectively, provided there are no selection ef- 
fects. 

The number of close frequency pairs is much higher than 
in the previous case and higher than the number of pairs ob- 
served in FG Vir. However, Fig. [7] is based on 160 frequen- 
cies, whereas there are 67 independent frequencies detected 
in FG Vir. Applying Eq. (|19[1 to FG Vir's frequency spec- 
trum yields 10.4 close frequency pairs, which is slightly over 
half the number of observed close frequ e ncy p airs in this 
star. Also, as pointed out in iLenz et al.1 <|2008T ). there arc 
7 frequency pairs where the separation is below 0.01 c/d. 
Equation (|19[l would yield 1.12 very close frequency pairs for 
FG Vir . Furthermore, according to iBreger fc Pamvatnvkbl 
(2006aj), close frequency pairs seem to cluster around ra- 
dial modes, whereas a number of pairs in Fig. [7] are close 
to the midpoints between consecutive radial modes. Hence, 
an astrophysical origin is still needed to account for this 
phenomena. 

An interesting feature appears in Fig. [7] Between ap- 
proximately 0.02 f^K and 0.05 Ok, few frequency pairs are 
detected. A simple explanation is that the rotation rate 
is sufficient to keep individual multiplet components far 
enough apart, but too small to cause frequency multiplets 
to overlap. Of course, including modes with higher £ values 
would probably reduce this gap by introducing new multi- 
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Figure 7. Upper panel: close frequency pairs in 2 Mq models as 
a function of the rotation rate for the frequency spectrum (n = 
1 — 10, I = — 3, m = — I to I). The horizontal dotted lines 
show the location of radial modes. Lower panel: number of close 
frequency pairs as a function of the rotation rate. The dashed line 
shows the number of close frequency pairs that is expected for a 
Poisson distribution. 



plets int o the frequency s pectr um. Ho wever, as was pointed 
out b y ILignieres et al.l (|2006l ) and ILignieres fc Georgeotl 
(2009), cancellation effects for such modes are more effec- 
tive at lower rotation rates, i.e. where close frequencies pairs 
are lacking. Table [2] gives a list of 5 Scuti stars where close 
frequencies have been reported. Interestingly, none of the 
stars fall in this gap. Of course, a larger number of stars 
would need to be analysed to see whether this gap remains 
or whether it is simply due to poor statistics. Furthermore, 
some of the sta rs in Table [21 such as BI CMi (IBreger et al.l 
l2002h . BV Cir (jMantegazza^ Poretti fc Zerbil 120011 ) and 44 
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Star HD v ■ sin i Reference 





number 


(km.s 1 ) 




BI CMi 


66853 


76 ±1 


Breeer et al. (20021 


XX Pyx 




52 ±2 


Handler et al. ( 19971 


V509 Per 


18878 


134 


Bush & Hintz (20081 


4 CVn 


107904 


112 ±3 a 


Bush & Hintz (20081 


FG Vir 


106384 


21.6 ± 0.3 b 


Zima et al. (20061 


BV Cir 


132209 


96.5 ± 1 


Mantcgazza ct al. (2001) 


BW Cnc 


73798 


200 ± 11 


Fossati et al. (20081 


44 Tau 


26322 


2±l c 


Zima et al. (20071 



thereby favouring an astrophysical origin to the phenom- 
ena. At lower rotation rates, where frequency multiplets do 
not overlap, a lack of close frequency pairs is observed in 
the calculated spectra, except at the lowest rotation rates 
where the rotational shift is small. Delta Scuti stars cur- 
rently reported as having close frequency pairs do not lie in 
this interval, but a more systematic study is required to see 
to what extent this is due to poor statistics, and to see what 
effects including gravity type modes has for more evolved 
stars. 



Table 2. Projected equatorial velocities of 5 Scuti stars 
with observed close fr eque ncy pairs (this l ist is based on 
iBreger fe Bischoj|2002l and lLenz et~aUl2008l V The first two 
columns identify the star. The third and fourth columns give 
the projected equatorial v elocity and the refe rence where 
this velocity is obtained. a iBush fe Hintzl J2008fl give a list of 
several values. The value given here is the average plus or mi- 
nus the standard deviation. b Alt hough the projected equatorial 
velocity is small, Zi ma et al. I te006l l estimate that the t rue equa- 
torial velocity is 66 ± 16 km.s -1 . c IZima et all d2007T) estimate 
that the true equatorial velocity is 3 ± 2 km.s -1 . 

Tau (|Zima et al.ll2007T ). are evolved. Consequently, their pul- 
sation spectra are likely to contain acoustic, gravity (or 
gravito-inertial) and possibly mixed modes. 
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6 CONCLUSION 

In this study, the effects of stellar rotation on acoustic pulsa- 
tion frequencies have been investigated using a second order 
perturbative and a two-dimensional approach. A compar- 
ison of the two shows that perturbative methods are valid 
for equatorial velocities up to 30 and 50 km.s -1 for a CoRoT 
long and short run, respectively. The associated validity do- 
mains close l y ma tch the domains previously obtained in 
iReese et al] f2006r i for polytropic models. These results show 
that perturbative methods, which are simpler to work with 
and less demanding numerically, will remain useful for a 
number of stars. Nonetheless, various qualitative and quan- 
titative differences appear between this approach and two- 
dimensional calculations at sufficient rotation rates. Among 
these is the pairing up of adjacent mo des in frequency mult i- 
plets. This phenomena, first noted in lEspinosa et all (|2004t ). 
also shows up in the two-dimensional calculations presented 
here but not in the second order perturbative calculations, in 
which the spacing between consecutive modes varies mono- 
tonically. The cause of this pairing up remains a mystery, 
given that different modes in a multiplet are not coupled 
due to their differing azimuthal orders. 

The case of observed close frequency pairs is then dis- 
cussed. In particular, it is shown that adjacent mode pairs in 
frequency multiplets do not come close enough together to 
provide a likely explanation, as oppos ed to what was previ- 
ously suggested (|Espinosa et al.|[2004l l. A systematic search 
for all of the close frequencies pairs in the calculated spectra 
was then carried out. Results showed that at sufficient ro- 
tation rates, the number of close pairs of low degree acous- 
tic modes matched what is expected based on a Poisson 
distribution. Even then, this does not match the number 
and distribution of close pairs observed in stars like FG Vir 



References 

Breger M., Bischof K. M., 2002, A&A, 385, 537 

Breger M. et al. 2002, MNRAS, 329, 531 

Breger M., Pamyatnykh A. A., 2006a, MNRAS, 368, 571 

Breger M., Pamyatnykh A. A., 2006b, Memorie della Soci- 

eta Astronomica Italiana, 77, 295 
Bush T. C, Hintz E. C, 2008, Astr. J., 136, 1061 
Christensen-Dalsgaard J., 2008a, ApSS, 316, 113 
Christensen-Dalsgaard J., 2008b, ApSS, 316, 13 
Clement M. J., 1981, ApJ, 249, 746 
Dziembowski W. A., Goode P. R., 1992, ApJ, 394, 670 
Eggleton P. P., Faulkner J., Flannery B. P., 1973, A&A, 

23, 325 

Espinosa F., Perez Hernandez F., Roca Cortes T., 2004, 
in ESA SP-559: SOHO 14 Helio- and Asteroseismology: 
Towards a Golden Future Oscillation Modes in Axially 
Symmetric Stars, pp 424-427 

Fossati L., Bagnulo S., Landstreet J., Wade G., Kochukhov 
O., Monier R., Weiss W., Gebran M., 2008, A&A, 483, 891 

Gough D. O., Thompson M. J., 1990, MNRAS, 242, 25 

Handler G. et al. 1997, MNRAS, 286, 303 

Iglesias C. A., Rogers F. J., 1996, ApJ, 464, 943 

Karami K., Christensen-Dalsgaard J., Pijpers F. P., Goupil 
M.-J., Dziembowski W. A., 2005, astro-ph/0502194 

Kjeldsen H., Christensen-Dalsgaard J., Handberg R., 
Brown T. M., Gilliland R. L., Borucki W. J., Koch D., 
2010, astro-ph.SR/1007.1816 

Kurucz R. L., 1991, in NATO ASIC Proc. 341: Stellar At- 
mospheres - Beyond Classical Models New opacity calcu- 
lations, p. 441 

Lenz P., Pamyatnykh A. A., Breger M., 2008, Journal of 

Physics Conference Series, 118, 012063 
Lignieres F., Georgeot B., 2008, Phys. Rev. E, 78, 016215 
Lignieres F., Georgeot B., 2009, A&A, 500, 1173 



On the effects of rotation on acoustic stellar pulsations 



Lignieres F., Rieutord M., Reese D., 2006, A&A, 455, 607 
Lovekin C. C, Deupree R. G., 2008, ApJ, 679, 1499 
MacGregor K. B., Jackson S., Skumanich A., Metcalfe 

T. S., 2007, ApJ, 663, 560 
Mantegazza L., Poretti E., Zerbi F. M., 2001, A&A, 366, 

547 

Michel E. et al. 2008, Communications in Asteroseismology, 
156, 73 

Ouazzani R.-M., Goupil M.-J., Dupret M.-A., Reese D., 
2009, Communications in Asteroseismology, 158, 283 

Reese D., Lignieres F., Rieutord M., 2006, A&A, 455, 621 

Reese D. R., MacGregor K. B., Jackson S., Skumanich A., 
Metcalfe T. S., 2009, A&A, 506, 189 

Saio H., 1981, ApJ, 244, 299 

Soufi F., Goupil M.-J., Dziembowski W. A., 1998, A&A, 
334, 911 

Suarez J. C, Goupil M. J., Reese D. R., Samadi R., 
Lignieres F., Rieutord M., Lochard J., 2010, ApJ, 721, 
537 

Yoshida S., Eriguchi Y., 2001, MNRAS, 322, 389 

Zima W., Lehmann H., Stiitz C, Ilyin I. V., Breger M., 

2007, A&A, 471, 237 
Zima W. et al. 2006, A&A, 455, 235 



